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ABSTRACT Mechanisms for electron transfer within microbial aggregates derived from an upflow anaerobic sludge blanket reac- 
tor converting brewery waste to methane were investigated in order to better understand the function of methanogenic consor- 
tia. The aggregates were electrically conductive, with conductivities 3-fold higher than the conductivities previously reported for 
dual-species aggregates of Geobacter species in which the two species appeared to exchange electrons via interspecies electron 
transfer. The temperature dependence response of the aggregate conductance was characteristic of the organic metallic-like con- 
ductance previously described for the conductive pili of Geobacter sulfurreducens and was inconsistent with electron conduction 
through minerals. Studies in which aggregates were incubated with high concentrations of potential electron donors demon- 
strated that the aggregates had no significant capacity for conversion of hydrogen to methane. The aggregates converted formate 
to methane but at rates too low to account for the rates at which that the aggregates syntrophically metabolized ethanol, an im- 
portant component of the reactor influent. Geobacter species comprised 25% of 16S rRNA gene sequences recovered from the 
aggregates, suggesting that Geobacter species may have contributed to some but probably not all of the aggregate conductivity. 
Microorganisms most closely related to the acetate-utilizing Methanosaeta concilii accounted for more than 90% of the se- 
quences that could be assigned to methane producers, consistent with the poor capacity for hydrogen and formate utilization. 
These results demonstrate for the first time that methanogenic wastewater aggregates can be electrically conductive and suggest 
that direct interspecies electron transfer could be an important mechanism for electron exchange in some methanogenic sys- 
tems. 

IMPORTANCE The conversion of waste organic matter to methane is an important bioenergy strategy, and a similar microbial 
metabolism of complex organic matter in anaerobic soils and sediments plays an important role in the global carbon cycle. Stud- 
ies with laboratory cultures have demonstrated that hydrogen or formate can serve as an electron shuttle between the microor- 
ganisms degrading organic compounds and methanogens. However, the importance of hydrogen and formate as intermediates 
in the conversion of organic matter to methane in natural communities is less clear. The possibility that microorganisms within 
some natural methanogenic aggregates may directly exchange electrons, rather than producing hydrogen or formate as an inter- 
mediary electron carrier, is a significant paradigm shift with implications for the modeling and design of anaerobic wastewater 
reactors and for understanding how methanogenic communities will respond to environmental perturbations. 
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Effective interspecies electron exchange is essential in methano- 
genic ecosystems (1-3). In order for complex organic matter 
to be converted to methane and carbon dioxide, the microbial 
community that metabolizes multicarbon organic compounds 
other than acetate requires an electron sink. Methanogenic micro- 
organisms can function as this sink, consuming electrons in the 
reduction of carbon dioxide to methane. Understanding the 
mechanisms of this electron exchange is key to modeling and/or 
manipulating methane production in natural methanogenic en- 
vironments, such as anaerobic soils and sediments, the gastroin- 



testinal systems of diverse animals, and anaerobic wastewater 
treatment systems. 

The first mechanism described for electron exchange in metha- 
nogenic systems was interspecies hydrogen transfer, in which mi- 
croorganisms that require an electron sink reduce protons to pro- 
duce hydrogen and the methanogens utilize hydrogen as an 
electron donor (1, 2, 4, 5). Formate may also act as an electron 
carrier (1, 5-7). 

A potential alternative is direct electron transfer between cells. 
Connections between cells via electrically conductive filaments 
were suggested as a mechanism for microbes to exchange elec- 
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FIG 1 Scanning electron micrographs of an entire aggregate (A) and higher magnification of the aggregate 
surface (B). 



trons (8-10). However, only one specific example of direct elec- 
tron transfer to a methanogen was proposed (9), and subsequent 
studies demonstrated that the filament connecting the fermenta- 
tive microorganism and the methanogen was a fiagellum, which is 
not expected to be conductive (2, 11). Therefore, direct electron 
transfer to the methanogen was unlikely (2). 

However, studies in which two Geobacter species were cocul- 
tured under conditions that required interspecies electron ex- 
change indicated that direct electron transfer is feasible (12). Ad- 
aptation of Geobacter metallireducens and Geobacter 
sulfurreducens to cooperatively metabolize ethanol resulted in the 
formation of large (1 to 2 mm in diameter) aggregates. Analysis of 
a mutation selected over the course of the laboratory evolution of 
the coculture, as well as studies in which specific mutations were 
introduced into G. sulfurreducens, suggested that the cells within 
the aggregates were directly exchanging electrons via a mechanism 
that involved the multiheme c-type cytochrome OmcS (12). Pre- 
vious studies (13) demonstrated that OmcS aligns along the con- 
ductive pili of G. sulfurreducens, and a similar localization was 
observed in the aggregates, suggesting that pili were also involved 
in the electron exchange. Furthermore, the aggregates were elec- 
trically conductive, demonstrating that electrons that G. metalli- 
reducens released into the extracellular matrix might readily flow 
to G. sulfurreducens ( 12). Studies of the capacity for hydrogen and 
formate metabolism suggested that they were not important elec- 
tron carriers between the two Geobacter species. 

The aggregates that the Geobacter species formed had a mor- 
phology similar to that of the aggregates that are found in upflow 
methanogenic wastewater digesters (14). Methanogenic digester 
aggregates are comprised of a diverse community of hydrolytic- 
fermentative bacteria, hydrogen-producing acetogenic bacteria, 
and methanogens, which cooperate to degrade complex organic 
compounds to methane and carbon dioxide (1, 2, 15, 16). One 
factor favoring aggregation is that it prevents cells from being 
washed out of the system (1, 14, 17-19). Furthermore, aggregation 
could enhance the exchange of hydrogen and formate between 
syntrophic partners (20, 21). 

An additional potential benefit of aggregation is that it may 
make it feasible for cells to directly exchange electrons (12, 22). 
Here we report on conductive properties, community composi- 
tion, and metabolism of wastewater aggregates consistent with 
direct exchange of electrons within the aggregates. 



RESULTS AND DISCUSSION 

General characteristics of the ag- 
gregates. The aggregates obtained 
from the industrial-scale reactor and 
the aggregates propagated in the lab- 
oratory reactor had similar appear- 
ances. They were black and roughly 
spherical (diameter, 0.5 to 2 mm), 
comprised primarily of tightly 
packed rods (Fig. 1). 

Aggregate conductance. The 
possibility that the aggregates might 
be electrically conductive, in a man- 
ner similar to that described for the 
aggregates formed by a coculture of 
Geobacter species (12), was investi- 
gated. Aggregates spanning a non- 
conductive gap between two gold 
electrodes exhibited a linear current-voltage response when a DC 
voltage was applied, consistent with ohmic conductivity (Fig. 2A). 
The conductivity of the aggregates from the industrial-scale reac- 
tor (6.1 ± 0.3 ju-S/cm [mean ± standard deviation]; n = 3) and 
that of the aggregates propagated in the laboratory reactor (7.2 ± 
3.0 /j,S/cm) were similar. They were more conductive than the 
previously described conductive aggregates comprised of a cocul- 
ture of two Geobacter species (1.4 ± 0.3 ju,S/cm). Conductive alu- 
minum beads treated in a similar manner had a high conductivity 
(11.4 ±0.1 /xS/cm), whereas artificial wastewater alone (Fig. 2A) 
or porous alginate beads in wastewater (Fig. 2B) did not exhibit an 
ohmic response. 

In order to further evaluate the nature of the conductance, 
aggregate conductance was measured as a function of temperature 
(Fig. 3). Upon cooling from room temperature, aggregate con- 
ductivity increased exponentially by more than an order of mag- 
nitude. An exponential increase in conductivity upon cooling is a 
characteristic signature of organic metals, such as polyaniline, 
whereas inorganic metals and minerals exhibit only linear in- 
creases in conductivity (23). Therefore, even though the aggre- 
gates contain inorganic iron that could conceivably contribute to 
their conductivity, the temperature dependence response rules 
this out as a major source of conductivity and suggests an organic 
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FIG 3 Temperature dependence of aggregate conductivity measured using a 
four-probe setup. 
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FIG 2 Current-voltage response of reactor aggregates (A) and aluminum and alginate control beads servin 
as positive and negative controls, respectively (B), in artificial wastewater. 



metal and hence a biological response. This temperature response 
of the aggregates is similar to that of conductive biofilms and 
purified pilus preparations of Geobacter species (24). The contin- 
ued exponential decrease in conductivity with further cooling 
(Fig. 3) can be attributed to intrinsic disorder present in the ag- 
gregate samples (23, 24). 

Experimental evaluation of electron sources for methane 
production. To further evaluate possible mechanisms for electron 
exchange within the aggregates, the potential for methane produc- 
tion from ethanol and potential intermediates of ethanol metab- 
olism were evaluated (Fig. 4). High concentrations of potential 
substrates were added in order to estimate the maximum rates at 
which the compounds might be metabolized. Methane produc- 
tion rates from ethanol exceeded those from other substrates 
(Fig. 4). There was substantial production of methane from ace- 
tate, consistent with the necessity of the community to metabolize 
the acetate expected to be produced from the ethanol and propi- 
onate in the reactor feed and the fact that acetate was directly 



provided in the reactor feed. The ad- 
dition of hydrogen did not stimulate 
methane production, indicating that 
the aggregates were not adapted to 
convert hydrogen to methane. 

The aggregates had some poten- 
tial to convert formate to methane. 
However, the combined maximum 
capacity for the aggregates to convert 
acetate and formate to methane was 
less than the capacity for conversion 
of ethanol to methane (Fig. 4). This 
was the case even though the concen- 
trations of acetate and formate pro- 
vided to the aggregates were much 
higher than those the methanogens 
were likely to experience during syn- 
trophic metabolism of the added 
ethanol. Thus, if formate was the 
dominant electron shuttle for etha- 
nol metabolism, combined rates of 
acetate and formate conversion to 
methane would be expected to be 
higher than the rates of methane 
production from added ethanol. The 
finding that the combined rates of 
methane production from added 
formate and acetate were, in fact, 
lower that from ethanol suggested 
that a mechanism other than inter- 
species hydrogen and formate trans- 
fer was required to account for the 
rapid metabolism of ethanol ob- 
served in the aggregates. 

Microbial community composi- 
tion. The 16S rRNA gene sequences 
in the aggregates were evaluated in 
order to gain further insight into ag- 
gregate metabolism. Rarefaction 
curves for the 16S rRNA gene librar- 
ies constructed with universal prim- 
ers or primers specific to the archaea 
reached an asymptote, suggesting adequate coverage for assessing 
the dominant members of the microbial community (Fig. 5). The 
microbial community of the aggregates propagated in the labora- 
tory reactor (Table 1 ) was similar to that of the aggregates directly 
harvested from the industrial-scale reactor (data not shown). The 
microbial community was phylogenetically diverse, as expected 
from previous studies (25-27). Bacteria accounted for 79% of the 
16S rRNA gene sequences recovered and included members of the 
Proteobacteria (33.9%), the Chloroflexi (15.8%), the Firmicutes 
(6.4%), and the Synergistetes (4.4%) (Table 1). 

The most abundant sequence recovered from the aggregates 
belonged to Geobacter species most closely related to Geobacter 
daltonii (Table 1). The finding that Geobacter sequences ac- 
counted for 25% of the sequences recovered is consistent with the 
finding that, in the absence of reactor disturbance, 20 to 30% of 
the 16S rRNA gene sequences recovered from aggregates treating 
brewery wastewater were in the family Desulfuromonadales, which 
includes Geobacter species (28). Although the most intensively 
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FIG 4 Methane production from the aggregates with different potential 
electron donors. Controls received no electron donor additions. Results are 
the means and standard deviations of data for six replicates for each treatment. 



studied physiological characteristic of Geobacter species is their 
capacity for extracellular electron transfer to metals and electrodes 
(22, 29, 30), they are also capable of syntrophic growth (12, 31- 
33). Syntrophic oxidation of ethanol and organic acids is their 
most likely role in brewery digester aggregates. 

Also abundant were sequences most closely related to Metha- 
nosaeta concilii (Table 1). The importance of this organism was 
further confirmed with the archaeal domain-specific primer set 
(Table 1). This is consistent with the known role of Methanosaeta 
species in initiating aggregate formation in digesters and their 
usual role as the primary methanogens converting acetate to 
methane in these systems (34, 35). 

Sequences that could be assigned to methanogens capable of 
utilizing hydrogen or formate were much less abundant than 
those of the acetate utilizers (Table 1), consistent with the low 
potential of the aggregates to convert these substrates to methane. 
The only 16S rRNA gene sequences recovered that were closely 
related to hydrogen/formate-utilizing methanogens were those 
related to Methanobacterium petrolearium, a hydrogen-utilizing 
methanogen (36). These sequences accounted for only 1.5% of the 
sequences recovered with the primers designed to recover se- 
quences from all microorganisms. The studies with both primer 
sets indicated that Methanobacterium sequences were less than 
10% as abundant as Methanosaeta sequences. These results fur- 
ther suggest that that neither hydrogen nor formate was an impor- 
tant electron shuttle for methane production. 



Implications. These results suggest that direct electron trans- 
fer could be an important mechanism for electron exchange 
within methanogenic wastewater treatment aggregates. The ag- 
gregates were even more conductive than the previously described 
aggregates of G. metallireducens and G. sulfurreducens, which mul- 
tiple lines of evidence suggested exchanged electrons via direct 
cell-to-cell electron transfer (12). Geobacter species were present 
in high abundance in the aggregates and could be expected to 
facilitate long-range electron transfer via their conductive pili (8, 
10, 37). However, the higher conductivity of the methanogenic 
aggregates suggests that aggregate constituents other than Geobac- 
ter species could have been contributing to aggregate conductiv- 
ity. 

Methane production via direct interspecies electron transfer, 
rather than via interspecies transfer of hydrogen or formate, 
would require that methanogens have the capacity to directly ac- 
cept electrons as an electron donor for methane production. The 
apparent ability of a Methanobacterium-like isolate (38) and a 
strain of Methanococcus maripaludis (39) to accept electrons di- 
rectly from metallic iron suggests that this is possible. Therefore, it 
seems likely that the Methanobacterium species in the methano- 
genic aggregates could have the capacity for directly accepting 
electrons from the conductive aggregate matrix. It has also been 
proposed that the archaea within anaerobic methane-oxidizing 
aggregates, which are related to Methanosaeta (40, 41), are also 
capable direct accepting electrons (1, 42) because physiological 
and modeling studies have suggested that interspecies hydrogen 
and formate transfer is unlikely within those aggregates (5, 43- 
47). Although Methanosaeta has been considered to comprise ex- 
clusively acetate -utilizing methanogens based on physiological 
studies, the genome sequences of Methanosaeta species contain 
genes encoding the enzymes required for the carbon dioxide re- 
duction pathway (41). Therefore, it cannot be ruled out that the 
highly abundant Methanosaeta species might also have the ability 
to directly accept electrons within aggregates. 

The formation of large, visible aggregates maybe only the most 
extreme instance of cells forming electrical contacts. The possibil- 
ity of cell-to-cell electron transfer in smaller methanogenic asso- 
ciations and in aggregates carrying out other types of synthrophic 
processes warrants further investigation. 

MATERIALS AND METHODS 

Operation of UASB reactor. Microbial aggregates were obtained from an 
upflow anaerobic sludge blanket (USAB) treating brewery waste at an 
Anheuser-Busch facility in Baldwinsville, NY. The industrial UASB reac- 
tor had an average chemical oxygen demand of 5.5 g liter - 1 and an average 
hydraulic retention time of 12 h. Aggregates were transported under an- 
aerobic conditions at room temperature before propagation in the labo- 
ratory. Aggregates for DNA extraction were stored at — 80°C prior to 
processing. 

The aggregates were propagated in the laboratory in a mesophilic 
(37°C) laboratory-scale (0.9 liters) UASB reactor (Fig. 6). The reactor was 
fed an artificial brewery wastewater, modified from the previously de- 
scribed synthetic brewery wastewater (48). No attempt was made to re- 
move dissolved oxygen from the artificial wastewater prior to its introduc- 
tion into the reactor, which simulates industrial conditions. The artificial 
wastewater composition (per liter) was as follows: ethanol, 2.5 ml; sodium 
propionate, 0.72 g; sodium acetate trihydrate, 1.8 g; urea, 0.23 g; KH 2 P0 4 , 
0.11 g; K 2 HP0 4 , 0.17 g; Na 2 S0 4> 0.05 g; MgCl 2 ■ 6H 2 0, 0.1 g; 
CaCl 2 ■ 2H 2 0, 0.05 g; trace element solution, 10 ml; vitamin solution, 
10 ml. The composition of the trace element solution (in grams per liter) 
was as follows: MnS0 4 • H 2 0, 0.5; FeS0 4 ■ 7H 2 0, 0.1; NiCL, ■ 6H 2 0, 0.04; 
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FIG 5 Rarefaction curves demonstrating high diversity coverage values for the 16S rRNA gene sequences 
recovered with the two primer sets employed. 



CoCl 2 ■ 6H 2 0, 0.048; ZnCl 2 , 0.13; CuS0 4 • 5H 2 0, 0.01; 
A1K(S0 4 ) 2 ■ 12H 2 0, 0.01; H 3 B0 3 , 0.01; Na 2 Mo0 4 ■ 2H 2 0, 0.025. The 
composition of the vitamin solution (in grams per liter) was as follows: 
biotin, 0.002; pantothenic acid, 0.005; B-12, 0.0001;p-aminobenzoic acid, 



0.005; thioctic acid (alpha lipoic), 0.005; 
nicotinic acid, 0.005; thiamine, 0.005; ri- 
boflavin, 0.005; pyridoxine HC1, 0.01; fo- 
lic acid, 0.002. The main carbon sources 
of the artificial wastewater (ethanol, pro- 
pionate, and acetate) created an approx- 
imate chemical oxygen demand of 
5.8 g liter -1 . The hydraulic retention 
time was 1 day. 

The reactor functioned under steady- 
state conditions throughout the course of 
the study. Volatile fatty acids were main- 
tained at low levels (0.18 ± 0.13 mM 
[mean ± standard deviation]), and 
chemical oxygen demand removal was 
96.7% ± 1.1%. Biogas production rates 
were 2,074 ± 341 milker -1 day -1 , witha 
methane composition of 78.9% ± 4.3%. 

Scanning electron microscopy. 
Scanning electron micrographs were 
taken of whole aggregates fixed with 2.5% 
glutaraldehyde for 4 h at 4°C, dehydrated 
and dried with hexamethyldisilazane 
(Sigma Aldrich) as previously described 
(49), and micrographs were taken as pre- 
viously described (12). 

Aggregate conductance. Two-probe 
electrical conductance measurements 
were performed with two gold electrodes 
separated by a 50-ixm nonconductive gap, as previously described (12). 
Aluminum and alginate beads and wastewater aggregates were suspended 
in artificial brewery wastewater before being placed on the gold electrodes, 



TABLE 1 Sequence abundances in 338F/907R and archaeal clone libraries with closest representative in GenBank" 



Library Phylum Closest relative'' No. of clones Relative abundance (%) 

338F/907R primers Actinobacteria Bacterium Ellin 6510 (92) 2 1.0 

Bacteroidetes Prolixibacter bellariivorans (89) 9 4.4 

Chloroflexi Longilinea arvoryzae (90%) 16 7.8 

Anaerobic bacterium MO-CFX2 (89) 2 1.0 

Bacterium JN18_A7_F* (96) 13 6.4 

Chloroflexi bacterium BL-DC-9 (86) 2 1.0 

Deltaproteobacteria Geobacter daltonii c (97) 50 24.5 

Syntrophobacter sulfatireducens (95) 11 5.4 

Syntrophorhabdus aromaticivorans (94) 3 1.5 

Syntrophobacter fumaroxidans MPOB (92) 3 1.5 

Syntrophus aciditrophicus (95) 2 1.0 

Euryarchaeota Methanobacterium petrolearium (97) 3 1.5 

Methanosaeta concilii (96) 40 19.6 

Firmicutes Clostridials sp. SM4/1 (94) 2 1.0 

Oryza sativa Indica (97) 1 0.5 

Syntrophomonas sp. TB-6 (97) 10 4.9 

Spirochaetes Spirochaeta sp. isolate TM-3 (93) 5 2.5 

Spirochaetes bacterium SA-10 (93) 3 1.5 

Spirochaeta sp. isolate TM-3 (93) 1 0.5 

Synergistetes Synergistetes bacterium 7WAY-8-7 (98) 1 0.5 

Aminomonas paucivorans (89) 8 3.9 

Thermotogae Thermotogales bacterium MesGlAg4.2.16S.B (98) 1 0.5 

Unclassified bacteria 16 7.8 

Arch21F/Arch958R primers Euryarchaeota Methanosaeta concilii (99) 102 87.8 

Methanobacterium petrolearium (99) 7 6.1 

Crenarchaeota Unclassified Desulfurococcales 7 6.1 



" The aggregates from the UASB reactor were used for analyses. Representative sequences comprising more than 10% of the sequences of their respective library are highlighted in 
grey. 

b Percentages of nucleotide similarity are shown in parentheses. 
c Formerly strain Geobacter FRC-32. 
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FIG 6 Diagram of the laboratory-scale UASB reactor. 



such that materials were spanning the nonconductive gap. Voltage was 
applied across the gap using a Keithley 2400 source meter. Voltage was 
scanned from -0.3 Vto +0.3 Vin steps of 0.025V and then from +0.3 V 
to —0.3 V. For each measurement, current was measured 10 s after setting 
the voltage to allow the exponential decay of the transient ionic current in 
the gap and to measure steady-state electronic current (50). The Lab- 
VIEW data acquisition system (National Instruments, Texas) was used to 
collect the data. Igor Pro (WaveMetrics Inc., Oregon) was used for data 
fitting and analysis. 

Four-probe electrical conductance measurements were performed as 
previously described for studies on microbial biofilms and pili (24). Four 
gold electrodes were separated by 50-u,m nonconductive gaps. Electrodes 
were fabricated using standard photolithography processing. The total 
electrode size was 1.25 cm by 1 cm. Aggregates were gently crushed to 
form a homogeneous layer and were placed on the four electrodes. The 
gentle crushing allowed aggregates to form a confluent film on all four 
electrodes and kept them from moving apart after drying. A source meter 
(Keithley 2400) was used to apply a fixed current between the outer of the 
four electrodes and to measure the potential drop between the two inner 
electrodes (51), by measuring the voltage for each current every second 
over a period of 100 s, after reaching the steady state. An additional high- 
impedance voltmeter (Keithley 2000) was used to record the output volt- 
age of the current source to calculate conductance (Lange and Mirsky 
[51]). 

The temperature dependence of conductivity was measured as previ- 
ously described (24) with a physical property measurement system 
(PPMS-6000; Quantum Design, California). The temperature was varied 
from room temperature (300 K) down to 210 K in steps of 5 K, and 
four-probe conductance was measured at each temperature. All temper- 
ature dependence experiments were performed in vacuum. 

Metabolism of ethanol and potential intermediates. In order to eval- 
uate the relative potential of the aggregates to metabolize ethanol and 
potential intermediates in ethanol metabolism to methane, 1-ml aliquots 
of aggregates taken from the lower portion of the laboratory reactor at 
steady state were added to anaerobic pressure tubes under N 2 -C0 2 (80:20) 
containing 10 ml of artificial brewery wastewater with the ethanol, propi- 



onate, and acetate omitted. NH 4 C1 (29 mM) was used in place of urea as a 
nitrogen source. The pH of the medium was adjusted to 7 after autoclav- 
ing with 1 M KOH. Ethanol (20 mM), acetate (20 mM), formate (20 mM), 
or hydrogen (0.7 atm as part of an H 2 -C0 2 [80:20] mix) were added to six 
replicate tubes for each treatment under strict anaerobic conditions. 
Methane production at 37°C was monitored by gas chromatography as 
described below. 

Community composition. Aggregates were frozen at — 80°C and 
ground to a fine powder under liquid nitrogen with a mortar and pestle 
before processing with Bio 101 FastDNA soil kits (MP Biomedicals, Ohio) 
according to the manufacturer's instructions. DNA was quantified with a 
NanoDrop spectrophotometer (Thermo, Fisher Scientific, Delaware). 

Portions of 16S rRNA gene sequences were amplified with the univer- 
sal primer set 338F (5' ACTCCTACGGGAGGCAGC 3') (52) and 907R 
(5' CCGTCAATTCCTTTRAGTTT 3') (53, 54) or the archaeal domain- 
specific primer set Arch21F (5' TTCCGGTTGATCCYGCCGGA 3') and 
Arch958R (5' TTCCGGTTGATCCYGCCGGA 3' ) (55). DNA was ampli- 
fied with a minicycler PTC 200 (Ml Research) starting at 95°C for 5 min, 
followed by 35 cycles consisting of denaturation (45 s at 95°C), annealing 
(45 s at 50°C for 338F/907R and 55°C for Arch21F/958R), and extension 
(1 min at 72°C) and a final extension at 72°C for 10 min. PCR products 
were purified with a gel extraction kit (Qiagen), and clone libraries were 
constructed with a Topo p2.1 cloning kit (Invitrogen, California) accord- 
ing to the manufacturer's instructions. Plasmid inserts from each clone 
library were sequenced with the M13F and M13R primer pair at the Uni- 
versity of Massachusetts Sequencing Facility. 

The 16S rRNA sequences were trimmed and assembled with the soft- 
ware package Geneious, version 5.3 pro, and were compared to sequences 
deposited in the GenBank DNA database using the BLAST algorithm (56) . 
16S rRNA gene sequences were assigned to phylogenetically consistent 
(>97% similarity) higher-order bacterial taxonomy groups based on a 
naive Bayesian rRNA classifier using the RDP Classifier software program 
(57). Alignments were achieved with the ClustalX v. 1.83 program (58) 
and manually proofread and corrected when necessary with ProSeq v. 2.9 
software (59) before performing rarefaction analyses. A distance matrix 
based on the alignment was constructed using DNAdist from the PHYLIP 
software program, version 3.6, with default parameters (60, 61). These 
pairwise distances served as input to the DOTUR software program 
(distance-based operational taxonomic unit [OTU] and richness), ver- 
sion 1.53 (62), for clustering the sequences into OTUs of 97% defined 
sequence similarity. These clusters served to generate rarefaction curves 
and to make calculations of Chaol (63) estimator and diversity coverage. 
To estimate the diversity coverage and determine if the number of clones 
from each library was sufficient, homologous coverage of each DNA li- 
brary was calculated using the equation C x = 1 — (N x /n), where N x is the 
number of unique sequences (singletons) and n is the total number of 
sequences (64). 

Analytical methods. Chemical oxygen demand was determined using 
Hach's method 8000 (Hach DR/890 colorimeter procedures manual). In 
brief, samples were digested with potassium dichromate for 2 h at 150°C 
to form green chromic ion (Cr 3+ ), which was quantified with a calibrated 
preprogrammed colorimeter (Hach DR-890). Volatile fatty acids were 
measured with a high-performance liquid chromatograph, as previously 
described (65). The ethanol concentration was measured with a gas chro- 
matograph equipped with a headspace sampler and a flame ionization 
detector (Clarus 600; PerkinElmer Inc., California). Biogas volume was 
measured by the water displacement method (66). The methane content 
was determined using gas chromatography with a flame ionization detec- 
tor (67). 

Nucleotide sequence accession numbers. The sequences determined 
in this work have been submitted to the EMBL database and assigned 
accession no. FR823513 to FR823761 for the 16S rRNA bacterial library 
and FR823762 to FR823877 for the 16S rRNA archaeal library. 
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